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Normalized  frequencies  are  computed  for  a  rectangular,  isotropic  plate  resting  on  elastic  supports.  The 
normalized  frequencies  are  determined  using  eigensensitivity  analysis,  which  approximates  the 
eigenparameters  in  a  Mandarin  series,  yielding  an  approximate  closed-form  expression.  One  benefit  of 
the  approximate  closed-form  expression  is  its  computational  efficiency  and  yet  another  is  its 
application  of  re-analysis.  Accuracy  of  the  approximate  expression  is  assessed  by  comparing  results 
with  the  widely  used  Rayleigh-Ritz  method  using  orthogonal  polynomials  and  beam  shape  functions  in 
both  approaches.  Consideration  for  a  variety  of  edge  conditions  is  given  through  a  combination  of 
simply  supported,  clamped  and  free  boundary  conditions.  Results  indicate  that  the  accuracy  of  higher 
frequencies  computed  by  the  sensitivity  approach  is  highly  dependent  upon  choice  of  basis  function. 
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1.  Background 

Point  supported  plates  are  plates  that  have  prescribed 
displacements  at  a  number  of  discrete  locations  within  its 
domain.  Rectangular,  circular  as  well  as  elliptical  plates  represent 
geometries  of  interest.  Rigid  point  supports  have  a  prescribed 
displacement  of  zero  while  elastically  point  supported  have 
displacements  dependent  of  the  stiffness  at  the  support.  Electro¬ 
nic  circuit  boards,  solar  panels,  and  concrete  slabs  represent 
applications  that  can  be  effectively  modeled  as  point  supported 
plates.  A  brief  overview  of  recent  contributions  is  presented. 
Altekin  [1]  consider  both  bucking  and  free  vibration  of  elliptical 
plates  which  have  point  supports  along  the  symmetric  diagonals. 
A  Ritz  approach  is  used  to  solve  for  the  fundamental  frequency 
and  critical  buckling  load.  Huang  et  al.  [2]  studied  the  free 
vibration  of  tapered  isotropic  plates  on  rigid  point  supports.  A 
Mindlin  plate  theory  is  used  and  a  Green’s  function  approach  is 
used  to  generate  the  governing  equations.  Yu  [3]  used  the  method 
of  superposition  established  by  Gorman  [4]  to  analyze  the  free 
vibration  of  cantilever  plates  containing  an  attached  mass.  Zhou 
et  al.  [5]  utilized  a  three-dimensional  theory  to  study  the 
frequency  analysis  of  both  isotropic  and  composites  plates.  A 
finite  layer  formulation  is  used  to  model  the  structure  and  a 
hybrid  basis  function  is  introduced  to  adequately  satisfy  the 
displacement  constraints  at  the  point  supports. 

Sensitivity  analysis  seeks  to  assess  the  effect  of  a  parameter  on 
the  response  of  a  system.  Application  of  sensitivity  analysis  covers 
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all  fields  of  economics,  business,  science,  mathematics,  and 
engineering.  Structural  eigensensitivity  analysis  provides  a  direct 
method  to  assess  the  effect  of  system  parameters  on  the 
eigenvalues,  typically  the  frequencies  and  buckling  loads.  An 
early  contributor  in  this  area  of  research  was  Hearmon  [6]  who 
developed  a  one-term  formula  to  approximate  the  fundamental 
frequency  of  orthotropic  plates.  Bert  [7]  studied  the  optimal 
design  of  composite  plates  for  maximum  fundamental  frequency. 
Barton  and  Reiss  [8]  provided  approximate  closed-formed 
formula  for  uni-axial  and  bi-axial  buckling  of  symmetric  compo¬ 
site  plates.  Recently  Barton  [9][10]  has  applied  this  technique  to 
investigate  the  bucking  of  isotropic  plates  subject  to  combined  in¬ 
plane  loading  and  the  thermal  buckling  of  composite  plates  with 
clamped-free  boundary  conditions. 

In  this  paper,  eigensensitivity  analysis  is  used  to  determine  an 
approximate  closed-form  expression  which  is  used  to  compute 
frequencies  of  a  square  isotropic  plate  resting  on  four  elastic  point 
supports.  A  combination  of  boundary  conditions  is  considered 
including  simply  supported,  clamped  and  free.  The  article  is 
organized  into  three  sections  beginning  the  problem  formulation, 
an  overview  of  the  eigensensivity  approach,  and  the  results  and 
discussion  section. 


2.  Problem  formulation 

There  are  many  ways  to  formulate  the  equations  that  govern 
the  free  vibration  of  the  elastically  point  supported  plate 
including  the  Newtonian  mechanics  and  the  principle  of  work- 
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energy.  Here  the  principle  of  minimum  potential  energy  is 
preferred. 

The  total  potential  energy  77  for  the  system  containing  elastic 
supports  and  discrete  masses  is  given  as 


n  =  U-T  (1) 

Here  U  contains  the  strain  energy  of  the  plate  and  the 
supporting  springs  given  by 


+ 


/  d^W 

J 


2 

-2(1 -d) 


kpS(x-Xp)3(y-yp)w^ 


p  =  1 


d^W  d^W 


(X,  y)dxdy 


(2) 


where  D  is  the  flexural  stiffness,  kp  is  the  stiffness  of  the  pth 
spring,  V  is  Poisson’s  ratio  and  w(x,y)  is  the  transverse  displace¬ 
ment.  The  kinetic  energy  consists  of  energy  from  the  plate  and 
from  added  discrete  masses.  Here  T  is  given  by 

Tpiate  =  JJ  w2(x,  y)dxdy+  JJ  mqSix-Xci)S(y-yq)w^ix,  y)dxdy 


Above  rriq  is  the  qth  discrete  mass.  In  the  paper,  all  stiffnesses  are 
taken  to  be  equal  as  kp  and  all  masses  are  taken  to  be  equal  as  m. 
The  displacement  w(x,y)  may  be  expressed  as 

w(x,  y)=Y.Y.  0  ©  (4) 

i  =  lj  =  1 

where  4>i{xla)  and  6j{ylb)  are  shape  function  in  the  x-  and 
y-directions,  respectively.  The  shape  functions  are  selected  to 
satisfy  kinematic  and  static  boundary  conditions.  The  class  of 
functions  can  be  beam  shape  or  polynomials.  In  fact.  Smith  et  al. 
[11]  considered  the  buckling  of  isotropic  plates  under  shear 
loading  and  solved  the  corresponding  problem  using  the  Ray- 
leigh-Ritz  method.  These  authors  considered  a  number  of 
orthogonal  polynomials  as  basis  functions  including  Chebyshev 
type-1  and  type-2,  Legendre,  Hermite,  and  Laguerre. 

The  governing  differential  equation  comes  by  extremizing 
Eq.  (1)  with  respect  to  the  kinematic  variable  w(x,y).  Substituting 
Eq.  (4)  into  Eq.  (1)  and  minimizing  the  total  potential  energy  with 
respect  to  the  coefficients  Wy,  results  in  the  eigenvalue  problem  of 

^  ^ijmn)^mn=0  (5) 

m  n 

where 

—  ^4  i^imbjn  +  +  ^miTrij)  +  2(1  —D)R^Ei^6jjj  +  R'^Bjfj^Cljn)  +  l'^7,jmn 

Mijmn  =  bjn  +  mF  ymn  (6) 


and  Cymn  =  1  4>i(>(p)dj(yp)4>mO(p)Sr,(yp)  is  the  product  of  the 

basis  functions  and  is  evaluated  where  the  springs  and  masses  are 
located.  Several  parameters  are  used  in  the  above  definitions 
including  R  which  is  the  aspect  ratio  given  as  ajb  and  a  set  of 
matrices  called  boundary  condition  matrices.  These  boundary 
condition  matrices  are  deflned  by 

Apm  =  4^"m)  Clqn  =  (9" q,  9” n) 

Bpm  =  i4^p,  4^"m)  bqn  =  (9q,  9n) 

Cpm  =  (0^  p?  0m)  =  {9''q,  9n)  ^  ^ 

Epm  =  (0''p?  0''m)  6pm  =  i9fp,  9frn) 

where  the  prime  represents  the  spatial  derivative  and  ( • ,  • )  is 
used  to  represent  the  L2  inner  product  on  [0,1].  The  above 


equations  for  Kymn  and  Mymn  are  the  quite  general.  They  are 
independent  of  any  particular  set  of  basis  functions.  Therefore 
they  can  be  used  if  (/)p  and  9q  are  kinematically  admissible 
polynomials,  beam  shape  functions,  or  any  other  set  of  kinema¬ 
tically  admissible  functions. 


3.  Sensitivity  analysis 


A  complete  presentation  of  the  sensitivity  approach  used  in 
this  paper  can  be  referenced  in  Barton  and  Reiss  [8].  Eq.  (5), 
written  as 

[K]{CCy}  =  A.y[M]{(Xy}  (8) 

provides  the  basis  to  apply  the  sensitivity  approach.  Here  we 
deflne  [ocy]  as  the  (i  j)th  component  of  the  eigenvector  and  is  its 
corresponding  eigenvalue.  An  approximate  expression  for  the 
eigenvalue  can  be  determined  by  introducing  parameters  5i 
and  S2  into  Eq.  (8)  and  considering 

[k(S,)]{ay(SuS2)}  =  ly(SuS2)[M{S2)]{0Cy(SuS2)}  (9) 

where 


[I<(S,)]  =  [I<d]+S,[AI<] 
[M(S2)]  =  [Md]+S2[AM] 


Here  [Kd]  and  [Md]  are  a  diagonal  matrices  obtained  from  [K]  and 
[M]  by  deleting  all  off-diagonal  elements;  [A/C]  and  [AM]  are 
matrices  which  have  zeros  on  the  diagonal  and  contain  only  the 
off-diagonal  elements  of  [/C]  and  [M],  respectively.  The  parameters 
Si  and  S2  take  on  values  of  either  0  or  1.  If  both  Si  and  S2=0,  the 
solution  to  Eq.  (9)  becomes  the  ratio  of  the  diagonal  elements  of 
the  stiffness  matrix  [Kq]  and  matrix  [M^].  If  both  Si  and  S2  =  l,  then 
the  original  eigenvalue  problem,  Eq.  (9)  is  recovered.  The  desired 
eigenvalue  Imn  is  obtained  by  expanding  Imn  in  a  Maclaurin  series 
about  (Si,S2)=(0,0)  and  evaluating  at  (Si,S2)=(l,l).  Thus 

Aij  =  l,j(l ,  1 )  £  l,j(0, 0) + 3lij(0, 0)  +  i  0)  (10) 


The  expressions  appearing  on  the  right-hand  side  of  Eq.  (10)  were 
presented  in  [8]  and  are 


k(0,  0)=JL 


5A,j(0,  0)  =  0 

[KyyAMklij-MyyAKklij] ' 
kdkl^ijij -kjij^ldkl 


(11) 


ijij  k^il^j 


Substituting  Eq.  (11)  into  Eq.  (10)  provides  the  desired 
quadratic  approximate  closed-form  expression  of 

12" 

(12) 


]  _  1 


^mJ  Myy 


[KyyAMjdy-MyyAKklij]  " 


kdldq  ^ijij  kjijn  ^klkl 


4.  Discussion 

In  this  section,  the  general  form  of  the  approximate  closed- 
form  expression,  Eq.  (11),  is  presented  for  combinations  of  simply 
supported,  clamped  or  free  boundary  conditions  and  for  two  types 
of  basis  function.  Boundary  conditions  speciflc  results  are  then 
evaluated  based  upon  chosen  basis  functions. 

Eq.  (12)  is  specialized  for  a  plate  consisting  of  four  springs  of 
equal  stiffness  and  up  to  four  equal  masses  can  be  determined  by 
first  evaluating  the  stiffness  and  mass  matrix  elements  given 
through  Eq.  (6)  and  then  substituting  these  results  into  Eq.  (12). 
Evaluating  Eq.  (6)  requires  determining  the  diagonal  matrices  for 
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both  [K]  and  [M]  which  are 

h  D 

%  =  (A,biy  +  2uR2c,,c,j  +  2(l-.)R2£,,e^+R4B..a^+ 

(13) 

and  computing  the  off-diagonal  matrices  for  both  [K]  and  [M]  are 

AKijki  =  vF?{CikCji  +  CkiCjj)  +  2(^-v)R^Eikeji  +  — Fyki 

AMp,=  ^rp,  (14) 

Substituting  both  Eqs.  (13)  and  (14)  into  the  general  expres¬ 
sion  given  in  Eq.  (12)  yields 


(16) 


Simplification  to  Eq.  (16)  for  the  case  of  no  added  discrete  masses 
occurs  by  setting  m=0  resulting  in 


cofjQ'^ph 

^D“ 


knO^ 


A,,bjj^2I^CuCjj^R^Buajj^^r^^ 


<[BAr’-EE 


2R^QiC,j+  '^rMij 
Ai,i,b„+2R^Q,Ci,+I^Bkkan+ 


Ajjbjj  +  2R^CijCjj+R‘*Bjiajj+  ^  r,j,j  x 


(17) 


cofja^ph 
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Aii^jj  +  2 CiiCjj  +  2(1  —  v)R^Eiiejj + R'^Buajj  ^  F ijij  x  Babjj  F  j^F ijij 


^jj  ^ -r  ^ V  A  - -r  iv  uuujj  -r  — ^  i  ijij 

1 


/<na2 


EEi  1  \Aibjj  +  2vR^CiiCjj  +  2(l-v)R^Eiiejj+R^Biiaj^ 


x[^rw,j„]  X  X  U2(Q,9,  +  QiC,j)  +  2(l-.)R2£,,ej,+  ^^r®  I  x 


^ +  2(1  —  v)R^Eiiejj  +  R^BiiQjj  +  ^  ijij 

m 


^  IJIJ 


vR^  (CikCji  +  QjC/j)  +  2(1  -  v)R^Eikeji  +  T  ^aj 


knd^ 


Akk^ii  +  ^^f^^CkkQi  +  2(^-v)R^Ekkeii+R^BkkCiii  +  T } 


klkl 


m 


X  |finb;j-+  I 


]  ~  Afibjj  +  2 QjCjj  +  2(1  —  v)R^Eiiejj + R^BuQjj  +  ^  /" yy  x  ^  ^ w/d j  | 


(15) 


Eq.  (15)  is  valid  for  any  combination  of  simply  supported, 
clamped  or  free  boundaries.  The  formula  requires  identifying  an 
admissible  basis  and  evaluating  the  boundary  matrices  given  in 
Eq.  (7).  Selecting  orthogonal  polynomials  as  the  basis  requires 
numerical  evaluation  of  the  boundary  matrices  and,  as  a  result,  no 
simplification  occurs  in  the  form  of  Eq.  (15).  If  free  boundary 
conditions  are  excluded,  then  the  following  results  holds 
Epm=-Cpm  and  epm=-Cpm.  Eq.  (15)  then  becomes 


cof^F^ph 


fr,  u  ^ 

'  jvj  J 
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m  y 

ijij  J 


kryO^ 


EE  \Ai^jjF2^C,Cjj+R^B,ajj+^F^ 


k^il^j 
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1  Td  K  ^  ^  r  1 
X  j^A/ynJ  X  ^Biibjj+-Fijij  X 


D  'Jy 

Aiibjj + 2R2CuCjj + R'^Buajj  +  y-  ryij 
2R2Q,.c,j+^r® 


+  \  Biibjj+ 


kpQ^ 


2R  CkiCy  +  E^klij 


k„a^ 


Biibjj+’^ry^  - 


knd^ 


A,^jj  +  2R^CiiCjj  +  R%iajj+  -F^F, 


and  if  both  added  discrete  mass,  m  =  0,  and  elastic  support 
stiffness,  kp=0,  Eq.  (17)  reduces  to 

cof:a'^ph 


D 


—  ^Aiibjj  -\-2R  CiiCjj-\-R  BiiCijj^  x  [Bjjby] 

CXD  CXD 

-'Y^'^{[2R^CkiCjiif}  X  {[Akkbii  +  2R^CkkCii  +  R'^BkkCiii] 
k^il^j 

x[Biibjj]-[AiibjjF2R^CiiCjj+R^Biiajj]  x  [Bkkbu]}-^}  (18) 


Specific  results  can  now  be  evaluated  for  a  variety  of  support 
conditions  including  simply  support  on  all  sides,  simply  sup¬ 
ported  on  opposite  sides,  and  clamped  on  all  sides. 


5.  Simply  supported  plates  (S-S-S-S) 

The  first  plate  to  be  considered  is  one  containing  any  number 
of  discrete  masses  and  equal  elastic  supports  such  that  the  ith 
spring  coordinates  lie  within  0  <  X/  <  a  and  0  <  yi  <  b.  Eor  a  plate 
simply  supported  on  all  sides,  the  boundary  condition  in  the  x  and 
y-direction  is  given  by 

w  =  0  Mxx  =  0  X  =  0,  a 
w  =  0  Myy  =  0  y  =  0,b 

Selecting  orthonormal  polynomials  [11]  as  the  basis  functions 
requires  numerically  evaluating  the  boundary  matrices.  These 
results  are  recalled  when  needed  to  evaluate  Eq.  (16).  The 
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Gramm-Schmidt  process  is  used  to  construct  an  ortho-normal  set 
after  assuming  the  initial  element  to  be 

(^o(x)  =  x-2x^+x^  0o(y)=y-2y^+y‘* 

Choosing  beam  shape  functions  allows  for  an  analytical  evalua¬ 
tion  of  the  boundary  matrices  which  can  be  directly  substituted 
into  Eq.  (16).  To  this  end,  select  the  basis  function  as 

(19) 

Evaluating  the  boundary  matrices  given  by  Eq.  (7)  provides 


(20) 


and  substituting  these  into  Eq.  (16)  yields  the  desired  approx¬ 
imate  closed-form  expression  for  simply  supported  boundary 
conditions  on  all  sides. 
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then  the  boundary  matrices,  Eq.  (7)  become 


Am  = 

^im  —  ^im 
Cjm  =  —P' 71^  Sir 


Cljn  —  Sjxi 
^jn  =  ^jn 


and  Eq.  (16)  becomes 

D 


cofjQ'^ph 


fn^-27t^eR^Cjj+R^2;+ 


■  1  +  — r 'I  '  — 


1 


(1  +(m/M)r,j,j)^ 


<EE1  1  \i^7t^-27t^i^R^Cjj+R^2f  + 


ksO^ 

D 


x[^r,„]  -2x 
rm 


(23) 


i*n*-2n^^R^Cjj+R^2]+ 


rm  ^  1  m  ^  1  ^n2i2  2c  KpO  r. 

X  jvj^wyJ  X  [l  +  jvj^yyJ  x  -2R  k  n  5/dCy+-^rfe„j 


1  +  — r-'T 


-2R't'it=4,ci+!^r,„ 


V>xi-2nVl?c,+K'>.U%-rm 


x['  +  S''®]- 

x[i  +  2r„,]} 


.4  t4  i<-psU 


r7i^-27iH^R^Cjj-TR^2j-T 
-1 


(24) 


A  simplification  to  Eq.  (21)  for  the  case  of  no  added  discrete 
masses  produces 

(D^:a^ph 


^  =7r^[i^  +  2i?2/2j2^^4j-4^ 


-4n^  X  EEr?M=f 


[(k^a^/n^D)rai^ 


[(/d-i4)  +  2i?2(fc2/2_i2j2)  +  ^?4(/4_j4)  +  (fc^a2/;,4D)(rfc,,,-r^^^^^^ 


Plates  which  are  simply  supported  in  the  x-direction  and  free  in  y- 
directions  have  boundary  conditions  given  by 

w  =  0  Mxx  =  0  X  =  0,  a 

^ — h  ^  o  o — h  (2 — r>)  ^  y  =  0,  b 

dy^  dx^  dy^  ox^dy 

and  following  beam  shape  functions  of 
/m7ix\ 


(22)  ^„g)=V2Sin(i^) 


and  if  both  added  discrete  mass  m=0  and  elastic  support  stiffness 
kp=0,  Eq.  (22)  reduces  to  the  well  known  equation  for  the 
frequency  of  a  simply  supported,  isotropic  plate. 

(ola'^ph  ..  „  „  „  .  . 

^  =71^  [i^  +  2RH^j^+R'^f] 


6.  Opposite  edges  simply  supported 

In  a  similar  manner,  an  expression  for  a  plate  with  two 
opposite  edges  simply  supported  in  the  x-direction  and  clamped 
in  the  y-direction  requires  identifying  the  following  beam  shape 
functions 


are  used  to  satisfy  them.  As  in  the  case  of  a  plate  with  clamped 
conditions  along  the  y-direction,  the  coefficients  defined  above 
hold  for  this  case  of  free  boundary  in  along  the  y-direction.  The 
boundary  matrices  give  through  Eq.  (8)  become 


2ny' 


2ny 


2ny 


Aim  = 

Cljn  —  Sjfi 

Rim  —  ^im 

^jn  =  ^jn 

Cim  = 

Cjn  =  (O'),  On) 

Eim  =  irnn^dim 

Cjn  —  i0^j->  Ofn) 

As  before  substituting  the  boundary  matrices  into  Eq.  (16) 
provides 


co?a4ph 


i'^n'^-2vR^fn^(ejj-hCjj)-h2R^i^n^ejj-hR^lf-h 


Please  cite  this  article  as:  Joseph  Watkins  R,  Barton  0  Jr..  Characterizing  the  vibration  of  an  elastically  point  supported  rectangular 
plate  using  eigensensitivity  analysis.  Thin  Walled  Struct  (2010),  doi:  10.1 01 6/j.tws.2009.1 1.005 


ARTICLE  IN  PRESS 


R.  Joseph  Watkins,  0.  Barton  Jr.  /  Thin-Walled  Structures  l  fllllj  lll-lll 


5 


r,  m  _  ] 


(bjj-h(m/M)rijijf 
'E^^{^^\i^^*-2oRW(ejj+Cjj)+2RhVejj+R*Xf+'^r^^^^^ 

X  ^^7z'‘-2uR^p7z^(ejj+Ci,)+2RWeii+l^lf+'^rij,jj 
rmn]  X  [bM+'^ri,i\  X  [^4l]  +  [bM+'^ni,]\^  X 


rm 

4m 


X  [bj/+  ^ryy]  -  +  ^^yu] 


For  all  other  combination  of  boundary  conditions,  clamped  in  the 
x-direction  and  free  in  y-direction  for  instance,  Eq.  (16)  will  be 
used  since  no  simplification  occurs  after  evaluating  it  with  the 
boundary  matrices. 


7.  Results 

In  this  section,  numerical  results  are  presented  for  the  derived 
approximate  closed-form  results  and  compared  to  results  generated 
using  the  Rayleigh-Ritz  method  for  the  elastically  supported  plate. 
Under  consideration  is  a  square  plate  30  in  x  30  in  x  0.5  in  with 
Young’s  modulus  E  and  Poisson’s  ratio  v  of  10  x  10®  psi  and  0.315, 
respectively.  Each  of  the  fours  springs  has  the  same  stiffness  of 
1301b/in.  Choice  of  spring  placement  and  configuration  is  infinite. 
Three  possible  plate  configurations  are  shown  in  Fig.  1  with 
corresponding  spring  coordinates  given  in  Table  1.  Configuration  A 
has  springs  placed  along  the  main  dialogues,  B  has  springs  centered 
on  each  edge  and  C  has  springs  lumped  at  the  center  of  the  plate. 
Tabulated  results  are  presented  for  configuration  A. 

Accuracy  of  the  approximate  formula  is  compared  with  the 
Ritz  method.  Frequencies  are  first  computed  using  orthogonal 
polynomials  as  basis  function  and  then  compared  to  results  using 
beam  shape  functions.  For  convenience,  introduce  the  normalized 
frequency  parameter 

_  co.a^  [ph 
'  2n  \  D 

which  identifies  the  ith  normalized  Ritz  frequency  and  ky  is  used 
to  identify  (ij)th  normalized  frequency  computed  using  the 
quadratic  approximate  expression. 

Table  2  presents  convergence  results  for  /<i  and  ku  for  all 
boundary  conditions  using  orthogonal  polynomials  as  basis 
functions  and  Table  3  presents  the  same  data  using  beam  shape 


functions  instead  as  the  basis  functions.  Eq.  (16)  was  used  to 
generate  the  approximate  closed-form  results.  Configuration  A 
without  any  added  discrete  masses  was  used  to  generate  results 
in  all  the  tables.  In  general  the  closed-form  expression  accurately 
predicts  the  fundamental  frequency  for  all  boundary  conditions 
using  either  type  of  basis  function.  The  largest  percent  difference 
in  the  two  results  occurs  for  the  SFSF  boundary  condition. 
Although  negligible,  these  are  0.5%  and  0.73%  using  the 
orthogonal  polynomials  and  beam  shape  functions,  respectively. 

Tables  4-8  present  the  first  five  frequencies  computed  by  both 
methods.  Again  both  orthogonal  polynomials  and  beam  shape 


Table  1 

Coordinates  for  spring  placement  (inches). 


Spring 

Configuration  A 

Configuration  B 

Configuration  C 

X 

Y 

X 

Y 

X 

Y 

1 

7.5 

7.5 

15 

5 

15 

15 

2 

22.5 

7.5 

25 

15 

15 

15 

3 

22.5 

22.5 

15 

25 

15 

15 

4 

7.5 

22.5 

5 

15 

15 

15 

Table  2 

Convergence 

polynomials. 

results 

for  Ritz  and 

closed-form 

expression  using  orthogonal 

N 

SSSS 

CCCC 

SCSC 

SFSF 

CFCF 

Ritz 

Quad 

Ritz 

Quad 

Ritz 

Quad 

Ritz 

Quad 

Ritz 

Quad 

1 

3.149 

3.149 

5.730 

5.730 

4.609 

4.609 

_ 

_ 

_ 

_ 

3 

3.148 

3.148 

5.728 

5.728 

4.608 

4.608 

2.565 

2.565 

4.214 

4.214 

5 

3.147 

3.147 

5.727 

5.727 

4.608 

4.608 

2.564 

2.564 

4.213 

4.213 

7 

3.147 

3.147 

5.727 

5.727 

4.608 

4.608 

2.550 

2.553 

4.194 

4.196 

9 

3.147 

3.147 

5.727 

5.727 

4.608 

4.608 

2.550 

2.543 

4.190 

4.176 

11 

3.147 

3.147 

5.727 

5.727 

4.608 

4.608 

2.550 

2.535 

4.188 

4.189 

Table  3 

Convergence 

functions. 

results 

for  Ritz  and 

closed-form  expression  using  beam  shape 

N 

SSSS 

CCCC 

SCSC 

SFSF 

CFCF 

Ritz 

Quad 

Ritz 

Quad 

Ritz 

Quad 

Ritz 

Quad 

Ritz 

Quad 

1 

3.158 

3.158 

5.752 

5.752 

4.624 

4.624 

2.858 

2.858 

4.382 

4.382 

3 

3.158 

3.158 

5.736 

5.735 

4.617 

4.617 

2.779 

2.773 

4.323 

4.318 

5 

3.158 

3.158 

5.733 

5.732 

4.617 

4.617 

2.754 

2.757 

4.305 

4.306 

7 

3.158 

3.158 

5.733 

5.732 

4.616 

4.616 

2.742 

2.752 

4.297 

4.303 

9 

3.158 

3.158 

5.732 

5.732 

4.616 

4.616 

2.735 

2.750 

4.292 

4.301 

11 

3.158 

3.158 

5.732 

5.731 

4.616 

4.616 

2.729 

2.749 

4.288 

4.300 

(x,,y,) 

(x,.y.) 

(g) 

X® 

Elastic  Springs 

(g) 

•(g) 

Configuration  A 


Configuration  B 


Configuration  C 


Fig.  1.  Spring  locations. 
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Table  4 

Comparison  of  normalized  Ritz  frequencies  with  closed-formed  approximate 
frequencies  for  SSSS  boundary  condition  (Hertz). 


j 

Xj 

Xv 

X3i 

Poly 

Beam 

Poly 

Beam 

Poly 

Beam 

Poly 

Beam 

1 

3.147 

3.158 

3.147 

3.158 

7.254 

7.867 

9.423 

15.71 

2 

7.854 

7.867 

7.254 

7.867 

11.80 

12.58 

15.68 

20.42 

3 

7.859 

7.867 

9.424 

15.71 

15.68 

20.42 

21.54 

28.27 

4 

12.56 

12.58 

- 

26.70 

8.622 

31.42 

20.66 

39.27 

5 

15.71 

15.71 

- 

40.84 

- 

45.55 

- 

53.41 

Table  5 

Comparison  of  normalized  Ritz  frequency  with  quadratic  approximate  equation 
for  CCCC  boundary  condition,  N=ll,  /<=1612.5. 


j 

Xj 

Xy 

Xii 

hj 

Poly 

Beam 

Poly 

Beam 

Poly 

Beam 

Poly 

Beam 

1 

5.727 

5.732 

5.727 

5.731 

11.67 

11.85 

20.89 

20.99 

2 

11.68 

11.85 

11.67 

11.85 

17.22 

17.58 

26.21 

26.43 

3 

11.68 

11.85 

20.89 

20.99 

26.21 

26.43 

34.97 

34.94 

4 

17.22 

17.58 

32.92 

33.52 

38.23 

38.80 

46.95 

47.05 

5 

21.94 

20.94 

46.83 

49.20 

53.30 

54.41 

61.22 

62.50 

Table  6 

Comparison  of  normalized  Ritz  frequency  with  quadratic  approximate  equation 
for  SCSC  boundary  condition. 


J 

Xi 

Xy 

Poly 

Beam 

Poly 

Beam 

Poly 

Beam 

Poly 

Beam 

1 

4.608 

4.616 

4.608 

4.616 

8.168 

8.724 

10.20 

16.28 

2 

8.713 

8.722 

11.03 

11.213 

14.74 

15.26 

18.35 

22.52 

3 

11.03 

11.21 

20.43 

20.55 

24.39 

24.63 

29.09 

3\.ll 

4 

15.05 

15.26 

32.57 

33.17 

36.80 

37.32 

42.38 

44.47 

5 

16.26 

16.27 

46.53 

48.91 

51.11 

53.14 

57.48 

60.33 

Table  7 

Comparison  of  normalized  Ritz  frequency  with  quadratic  approximate  equation 
for  SFSF  boundary  condition  using  beam  shape  functions. 


j 

Xj 

hj 

1 

2.729 

2.749 

7.363 

15.17 

2 

7.352 

7.204 

11.48 

18.98 

3 

7.545 

10.25 

13.80 

20.67 

4 

12.06 

17.61 

20.32 

25.65 

5 

14.27 

22.20 

24.61 

29.63 

Table  8 

Comparison  of  normalized  Ritz  frequency  with  quadratic  approximate  equation 
for  CFCF  boundary  condition  using  beam  shape  functions. 


j 

Xj 

kij 

hj 

1 

4.29 

4.30 

10.81 

20.09 

2 

8.41 

8.06 

14.21 

23.30 

3 

10.80 

10.87 

16.16 

24.70 

4 

14.76 

17.94 

21.92 

28.88 

5 

14.84 

22.46 

25.97 

32.53 

functions  are  used  as  basis  functions.  All  results  are  generated  for 
configuration  A.  The  number  of  terms  taken  in  the  displacement 
expansion  corresponds  to  N=ll.  Frequencies  from  the  Ritz 
analysis  are  ordered  beginning  with  the  fundamental  frequency, 
the  lowest  numerical  result.  In  order  to  compare  frequencies 


computed  using  the  closed-form  expression,  several  modes 
corresponding  combinations  of  indices  i  and  j  must  be 
computed.  For  the  most  part,  /<i  =  kn,  /<2  =  /<i2.  ks  =  fei,  k4  =  (<22, 
and  ks  =  ksi. 

Table  4  provides  normalized  frequencies  for  the  simply 
supported  plate.  The  first  column  report  frequencies  computed 
using  the  Ritz  approach  using  both  types  of  basis  functions  while 
columns  two  through  four  report  frequencies  using  the  closed- 
form  expressions  given  in  Eq.  (16),  for  orthogonal  polynomials 
and  Eq.  (21)  for  beam  shape  functions.  Some  frequencies 
computed  using  the  closed-form  equation  with  orthogonal 
polynomials,  such  as  ki4,  are  not  real.  This  is  due  to  two  factors. 
First  (5^2i4,  which  contains  products  of  the  off-diagonal  elements 
of  [/C],  are  not  small  and  are  greater  than  the  zeroth-order  term. 
Secondly,  the  boundary  matrix  (p'UX  which  contributes  to 
the  off-diagonal  components  of  K,  is  not  small.  When  beam  shape 
functions  are  used  necessarily  0' j)  =  fifSy  where  jdi  is  the 
beam’s  ith  natural  frequency.  For  this  boundary  condition,  the 
approximate  closed-form  expression  Eq.  (21)  which  uses  beam 
shape  functions  compute  more  accurately  frequencies  when 
compared  with  the  Ritz  approach  using  beam  shape  functions 
than  does  the  closed-form  expression  Eq.  (16)  which  uses 
orthogonal  polynomials.  The  difference  in  results  is  more 
exaggerated  when  considering  higher  modes.  From  Table  4  using 
orthogonal  polynomials,  the  fourth  Ritz  mode,  /q,  is  12.56  Flz  and 
the  fourth  mode  k22  predicted  by  Eq.  (18)  is  11.80  Flz  resulting  in 
a  6.1%  difference.  Comparing  the  same  cases  using  beam  shape 
functions  shows  no  difference. 

Table  5  presents  results  for  a  plate  fully  clamped  on  all  sides. 
Eq.  (16)  was  used  to  generate  results  for  both  types  of  basis 
functions.  As  before,  the  closed-form  expression  better  predicts 
the  Ritz  frequencies  using  beam  shape  function  than  when  using 
orthogonal  polynomials.  A  percent  difference  of  4.7%  occurs  in  the 
comparison  of  the  fifth  frequency  using  polynomials  and  only 
0.24%  occurs  when  using  beam  shape  functions.  Table  6  shows  the 
first  of  three  results  for  mixed  boundary  conditions,  simply 
supported  in  the  x-direction  and  clamped  in  the  y-direction.  Eq. 
(23)  is  used  to  generate  all  numerical  results.  The  largest  percent 
difference  occurs  again  when  using  the  polynomials  as  basis.  Here 
a  12.9%  difference  occurs  in  the  computation  of  the  fifth 
frequency.  Percent  differences  of  6.3%  and  2.1%,  respectively 
occur  in  the  second  and  fourth  frequencies.  The  maximum 
percent  difference  of  0.06%  occurs  in  the  fifth  frequency  using 
beam  shape  functions. 

Tables  7  and  8  show  results  for  the  remaining  two  mixed 
boundary  condition  cases  which  includes  a  free  edge  in  the  y- 
direction.  Previous  results  indicate  that  beam  shape  functions  are 
better  suited  for  the  closed-form  expression,  especially  when 
considering  higher  frequencies.  Therefore  only  beam  shape  function 
results  are  shown.  For  the  SFSF  boundary  condition.  Table  7,  ki  is 
approximated  by  kn  and  differs  by  0.73%,  k2  is  approximate  by  ki2 
and  differs  by  2.01%,  ks  is  approximated  by  k2i  differing  by  2.4%,  k4  is 
approximated  by  k22  differing  by  4.8%,  and  ks  is  approximated  by  fes 
and  differs  by  3.3%.  For  the  SFSF  boundary  condition.  Table  7,  ki  is 
approximated  by  kn  and  differs  by  0.23%,  ks  is  approximate  by  ku 
and  differs  by  4.2%,  ks  is  approximated  by  ksi  differing  by  0.09%,  lu  is 
approximated  by  k22  differing  by  3.7%,  and  ks  is  approximated  by  kss 
and  differs  by  8.9%. 


8.  Conclusion 

In  this  paper  approximate  closed-form  formulas  were  devel¬ 
oped  to  determine  the  vibration  of  elastically  supported  plates.  A 
combination  of  edge  support  conditions  were  investigated 
indcluding  simply  supported,  clamped  and  free.  Both  orthogonal 
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polynomials  and  beam  shape  functions  were  used  as  admissible 
basis  functions.  When  compared  to  Rayleigh-Ritz  results,  the 
computed  fundamental  frequency,  for  all  boundary  conditions, 
using  the  approximate  closed-form  expressions  are  excellent.  The 
accuracy  of  higher  frequencies,  computed  using  the  approximated 
closed-form  expression  is  dependent  upon  boundary  condition 
and  choice  of  basis  functions.  A  maximum  percent  difference  of 
0.06%  occurs  using  beam  shape  functions  while  12.9%  occurs 
using  the  orthogonal  polynomials  for  the  SCSC  boundary  condi¬ 
tion.  For  other  mixed  boundary  conditions,  which  includes  a  free 
edge,  percent  differences  in  the  fifth  mode  of  4.8%  and  8.9%  occurs 
for  the  SFSF  and  CFCF  boundary  conditions,  respectively. 
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